A framework for conducting time-varying genome-wide association studies: An application to body mass index across childhood in six multiethnic cohorts.

Genetic effects on changes in human traits over time are understudied and may have important pathophysiological impact. We propose a framework that enables data quality control, implements mixed models to evaluate trajectories of change in traits, and estimates phenotypes to identify age-varying genetic effects in genome-wide association studies (GWASs). Using childhood body mass index (BMI) as an example, we included 71,336 participants from six cohorts and estimated the slope and area under the BMI curve within four time periods (infancy, early childhood, late childhood and adolescence) for each participant, in addition to the age and BMI at the adiposity peak and the adiposity rebound. GWAS on each of the estimated phenotypes identified 28 genome-wide significant variants at 13 loci across the 12 estimated phenotypes, one of which was novel (in DAOA) and had not been previously associated with childhood or adult BMI. Genetic studies of changes in human traits over time could uncover novel biological mechanisms influencing quantitative traits.

A tool to assess risk of bias in non-randomized follow-up studies of exposure effects (ROBINS-E).

BACKGROUND: Observational epidemiologic studies provide critical data for the evaluation of the potential effects of environmental, occupational and behavioural exposures on human health. Systematic reviews of these studies play a key role in informing policy and practice. Systematic reviews should incorporate assessments of the risk of bias in results of the included studies. OBJECTIVE: To develop a new tool, Risk Of Bias In Non-randomized Studies - of Exposures (ROBINS-E) to assess risk of bias in estimates from cohort studies of the causal effect of an exposure on an outcome. METHODS AND RESULTS: ROBINS-E was developed by a large group of researchers from diverse research and public health disciplines through a series of working groups, in-person meetings and pilot testing phases. The tool aims to assess the risk of bias in a specific result (exposure effect estimate) from an individual observational study that examines the effect of an exposure on an outcome. A series of preliminary considerations informs the core ROBINS-E assessment, including details of the result being assessed and the causal effect being estimated. The assessment addresses bias within seven domains, through a series of 'signalling questions'. Domain-level judgements about risk of bias are derived from the answers to these questions, then combined to produce an overall risk of bias judgement for the result, together with judgements about the direction of bias. CONCLUSION: ROBINS-E provides a standardized framework for examining potential biases in results from cohort studies. Future work will produce variants of the tool for other epidemiologic study designs (e.g. case-control studies). We believe that ROBINS-E represents an important development in the integration of exposure assessment, evidence synthesis and causal inference.

Using the global randomization test as a Mendelian randomization falsification test for the exclusion restriction assumption.

Mendelian randomization may give biased causal estimates if the instrument affects the outcome not solely via the exposure of interest (violating the exclusion restriction assumption). We demonstrate use of a global randomization test as a falsification test for the exclusion restriction assumption. Using simulations, we explored the statistical power of the randomization test to detect an association between a genetic instrument and a covariate set due to (a) selection bias or (b) horizontal pleiotropy, compared to three approaches examining associations with individual covariates: (i) Bonferroni correction for the number of covariates, (ii) correction for the effective number of independent covariates, and (iii) an r2 permutation-based approach. We conducted proof-of-principle analyses in UK Biobank, using CRP as the exposure and coronary heart disease (CHD) as the outcome. In simulations, power of the randomization test was higher than the other approaches for detecting selection bias when the correlation between the covariates was low (r2 < 0.1), and at least as powerful as the other approaches across all simulated horizontal pleiotropy scenarios. In our applied example, we found strong evidence of selection bias using all approaches (e.g., global randomization test p < 0.002). We identified 51 of the 58 CRP genetic variants as horizontally pleiotropic, and estimated effects of CRP on CHD attenuated somewhat to the null when excluding these from the genetic risk score (OR = 0.96 [95% CI: 0.92, 1.00] versus 0.97 [95% CI: 0.90, 1.05] per 1-unit higher log CRP levels). The global randomization test can be a useful addition to the MR researcher's toolkit.

Protein identification for stroke progression via Mendelian Randomization in Million Veteran Program and UK Biobank.

Background: Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke. Methods: We performed genome-wide association studies (GWAS) for subsequent major adverse cardiovascular events (MACE) (Ncases=51,929, Ncntrl=39,980) and subsequent arterial ischemic stroke (AIS) Ncases=45,120, Ncntrl=46,789) after first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (pQTLs) to determine the effect of 1,463 plasma protein abundances on subsequent MACE using Mendelian randomization (MR). Results: Two variants were significantly associated with subsequent cardiovascular events: rs76472767 (OR=0.75, 95% CI = 0.64-0.85, p= 3.69×10-08) with subsequent AIS and rs13294166 (OR=1.52, 95% CI = 1.37-1.67, p=3.77×10-08) with subsequent MACE. Using MR, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 (effect OR= 0.77, 95% CI = 0.66-0.88, adj. p=0.05), and TNFRSF14 (effect OR=1.42, 95% CI = 1.24-1.60, adj. p=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation. Conclusions: We found evidence that two proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.

Meta-regression of genome-wide association studies to estimate age-varying genetic effects.

Fixed-effect meta-analysis has been used to summarize genetic effects on a phenotype across multiple Genome-Wide Association Studies (GWAS) assuming a common underlying genetic effect. Genetic effects may vary with age (or other characteristics), and not allowing for this in a GWAS might lead to bias. Meta-regression models between study heterogeneity and allows effect modification of the genetic effects to be explored. The aim of this study was to explore the use of meta-analysis and meta-regression for estimating age-varying genetic effects on phenotypes. With simulations we compared the performance of meta-regression to fixed-effect and random -effects meta-analyses in estimating (i) main genetic effects and (ii) age-varying genetic effects (SNP by age interactions) from multiple GWAS studies under a range of scenarios. We applied meta-regression on publicly available summary data to estimate the main and age-varying genetic effects of the FTO SNP rs9939609 on Body Mass Index (BMI). Fixed-effect and random-effects meta-analyses accurately estimated genetic effects when these did not change with age. Meta-regression accurately estimated both main genetic effects and age-varying genetic effects. When the number of studies or the age-diversity between studies was low, meta-regression had limited power. In the applied example, each additional minor allele (A) of rs9939609 was inversely associated with BMI at ages 0 to 3, and positively associated at ages 5.5 to 13. Our findings challenge the assumption that genetic effects are consistent across all ages and provide a method for exploring this. GWAS consortia should be encouraged to use meta-regression to explore age-varying genetic effects.

Multiple imputation strategies for missing event times in a multi-state model analysis.

In clinical studies, multi-state model (MSM) analysis is often used to describe the sequence of events that patients experience, enabling better understanding of disease progression. A complicating factor in many MSM studies is that the exact event times may not be known. Motivated by a real dataset of patients who received stem cell transplants, we considered the setting in which some event times were exactly observed and some were missing. In our setting, there was little information about the time intervals in which the missing event times occurred and missingness depended on the event type, given the analysis model covariates. These additional challenges limited the usefulness of some missing data methods (maximum likelihood, complete case analysis, and inverse probability weighting). We show that multiple imputation (MI) of event times can perform well in this setting. MI is a flexible method that can be used with any complete data analysis model. Through an extensive simulation study, we show that MI by predictive mean matching (PMM), in which sampling is from a set of observed times without reliance on a specific parametric distribution, has little bias when event times are missing at random, conditional on the observed data. Applying PMM separately for each sub-group of patients with a different pathway through the MSM tends to further reduce bias and improve precision. We recommend MI using PMM methods when performing MSM analysis with Markov models and partially observed event times.

Investigation of genetic determinants of cognitive change in later life.

Cognitive decline is a major health concern and identification of genes that may serve as drug targets to slow decline is important to adequately support an aging population. Whilst genetic studies of cross-sectional cognition have been carried out, cognitive change is less well-understood. Here, using data from the TOMMORROW trial, we investigate genetic associations with cognitive change in a cognitively normal older cohort. We conducted a genome-wide association study of trajectories of repeated cognitive measures (using generalised estimating equation (GEE) modelling) and tested associations with polygenic risk scores (PRS) of potential risk factors. We identified two genetic variants associated with change in attention domain scores, rs534221751 (p = 1 × 10-8 with slope 1) and rs34743896 (p = 5 × 10-10 with slope 2), implicating NCAM2 and CRIPT/ATP6V1E2 genes, respectively. We also found evidence for the association between an education PRS and baseline cognition (at >65 years of age), particularly in the language domain. We demonstrate the feasibility of conducting GWAS of cognitive change using GEE modelling and our results suggest that there may be novel genetic associations for cognitive change that have not previously been associated with cross-sectional cognition. We also show the importance of the education PRS on cognition much later in life. These findings warrant further investigation and demonstrate the potential value of using trial data and trajectory modelling to identify genetic variants associated with cognitive change.

Two sample Mendelian Randomisation using an outcome from a multilevel model of disease progression.

Identifying factors that are causes of disease progression, especially in neurodegenerative diseases, is of considerable interest. Disease progression can be described as a trajectory of outcome over time-for example, a linear trajectory having both an intercept (severity at time zero) and a slope (rate of change). A technique for identifying causal relationships between one exposure and one outcome in observational data whilst avoiding bias due to confounding is two sample Mendelian Randomisation (2SMR). We consider a multivariate approach to 2SMR using a multilevel model for disease progression to estimate the causal effect an exposure has on the intercept and slope. We carry out a simulation study comparing a naïve univariate 2SMR approach to a multivariate 2SMR approach with one exposure that effects both the intercept and slope of an outcome that changes linearly with time since diagnosis. The simulation study results, across six different scenarios, for both approaches were similar with no evidence against a non-zero bias and appropriate coverage of the 95% confidence intervals (for intercept 93.4-96.2% and the slope 94.5-96.0%). The multivariate approach gives a better joint coverage of both the intercept and slope effects. We also apply our method to two Parkinson's cohorts to examine the effect body mass index has on disease progression. There was no strong evidence that BMI affects disease progression, however the confidence intervals for both intercept and slope were wide.

Modeling Risk Factors for Intraindividual Variability: A Mixed-Effects Beta-Binomial Model Applied to Cognitive Function in Older People in the English Longitudinal Study of Ageing.

Cognitive functioning in older age profoundly impacts quality of life and health. While most research on cognition in older age has focused on mean levels, intraindividual variability (IIV) around this may have risk factors and outcomes independent of the mean value. Investigating risk factors associated with IIV has typically involved deriving a summary statistic for each person from residual error around a fitted mean. However, this ignores uncertainty in the estimates, prohibits exploring associations with time-varying factors, and is biased by floor/ceiling effects. To address this, we propose a mixed-effects location scale beta-binomial model for estimating average probability and IIV in a word recall test in the English Longitudinal Study of Ageing. After adjusting for mean performance, an analysis of 9,873 individuals across 7 (mean = 3.4) waves (2002-2015) found IIV to be greater at older ages, with lower education, in females, with more difficulties in activities of daily living, in later birth cohorts, and when interviewers recorded issues potentially affecting test performance. Our study introduces a novel method for identifying groups with greater IIV in bounded discrete outcomes. Our findings have implications for daily functioning and care, and further work is needed to identify the impact for future health outcomes.

Young Adult ADHD Symptoms in the General Population and Neurocognitive Impairment.

OBJECTIVE: Neurocognitive impairments are associated with child and adult ADHD in clinical settings. However, it is unknown whether adult ADHD symptoms in the general population are associated with the same pattern of cognitive impairment. We examined this using a prospective, population-based cohort spanning birth to age 25 years. METHODS: We examined associations between self-reported adult ADHD symptoms and cognitive task performance (attention and response inhibition) in adulthood and childhood. RESULTS: Self-rated ADHD symptoms at age 25 were associated with poorer performance in age 25 cognitive tasks capturing ADHD-related functioning (attention B = -0.03, 95% CI [0.05, -0.01], p = .005; response inhibition B = -0.03, 95% CI [-0.05, -0.01], p = .002). CONCLUSIONS: Neurocognitive impairments linked to adult ADHD symptoms in the general population, are similar to those found in people with childhood ADHD symptoms and are consistent with findings in adult ADHD clinical samples.

An epigenome-wide analysis of DNA methylation, racialized and economic inequities, and air pollution.

Importance: DNA methylation (DNAm) provides a plausible mechanism by which adverse exposures become embodied and contribute to health inequities, due to its role in genome regulation and responsiveness to social and biophysical exposures tied to societal context. However, scant epigenome-wide association studies (EWAS) have included structural and lifecourse measures of exposure, especially in relation to structural discrimination. Objective: Our study tests the hypothesis that DNAm is a mechanism by which racial discrimination, economic adversity, and air pollution become biologically embodied. Design: A series of cross-sectional EWAS, conducted in My Body My Story (MBMS, biological specimens collected 2008-2010, DNAm assayed in 2021); and the Multi Ethnic Study of Atherosclerosis (MESA; biological specimens collected 2010-2012, DNAm assayed in 2012-2013); using new georeferenced social exposure data for both studies (generated in 2022). Setting: MBMS was recruited from four community health centers in Boston; MESA was recruited from four field sites in: Baltimore, MD; Forsyth County, NC; New York City, NY; and St. Paul, MN. Participants: Two population-based samples of US-born Black non-Hispanic (Black NH), white non-Hispanic (white NH), and Hispanic individuals (MBMS; n=224 Black NH and 69 white NH) and (MESA; n=229 Black NH, n=555 white NH and n=191 Hispanic). Exposures: Eight social exposures encompassing racial discrimination, economic adversity, and air pollution. Main outcome: Genome-wide changes in DNAm, as measured using the Illumina EPIC BeadChip (MBMS; using frozen blood spots) and Illumina 450k BeadChip (MESA; using purified monocytes). Our hypothesis was formulated after data collection. Results: We observed the strongest associations with traffic-related air pollution (measured via black carbon and nitrogen oxides exposure), with evidence from both studies suggesting that air pollution exposure may induce epigenetic changes related to inflammatory processes. We also found suggestive associations of DNAm variation with measures of structural racial discrimination (e.g., for Black NH participants, born in a Jim Crow state; adult exposure to racialized economic residential segregation) situated in genes with plausible links to effects on health. Conclusions and Relevance: Overall, this work suggests that DNAm is a biological mechanism through which structural racism and air pollution become embodied and may lead to health inequities.

Multiple imputation of missing data under missing at random: including a collider as an auxiliary variable in the imputation model can induce bias.

Epidemiological studies often have missing data, which are commonly handled by multiple imputation (MI). In MI, in addition to those required for the substantive analysis, imputation models often include other variables ("auxiliary variables"). Auxiliary variables that predict the partially observed variables can reduce the standard error (SE) of the MI estimator and, if they also predict the probability that data are missing, reduce bias due to data being missing not at random. However, guidance for choosing auxiliary variables is lacking. We examine the consequences of a poorly chosen auxiliary variable: if it shares a common cause with the partially observed variable and the probability that it is missing (i.e., it is a "collider"), its inclusion can induce bias in the MI estimator and may increase the SE. We quantify, both algebraically and by simulation, the magnitude of bias and SE when either the exposure or outcome is incomplete. When the substantive analysis outcome is partially observed, the bias can be substantial, relative to the magnitude of the exposure coefficient. In settings in which a complete records analysis is valid, the bias is smaller when the exposure is partially observed. However, bias can be larger if the outcome also causes missingness in the exposure. When using MI, it is important to examine, through a combination of data exploration and considering plausible casual diagrams and missingness mechanisms, whether potential auxiliary variables are colliders.

Childhood socioeconomic position and sex-specific trajectories of metabolic traits across early life: prospective cohort study.

BACKGROUND: Socioeconomic inequalities in cardiovascular disease risk begin early in life and are more pronounced in females than males later in life. Causal atherogenic traits explaining this are not well understood. We explored sex-specific associations between childhood socioeconomic position (SEP) and molecular measures of systemic metabolism across early life. METHODS: Data were from the Avon Longitudinal Study of Parents and Children (ALSPAC), a population-based birth cohort in southwest England. Pregnant women with an expected delivery date between 1991 and 1992 were invited to participate. Maternal education was the primary indicator of SEP. Concentrations of 148 metabolic traits from targeted metabolomics (nuclear magnetic resonance spectroscopy) from research clinics at ages 7, 15, 18 and 25 years were analysed. The sex-specific slope index of inequality (SII) in trajectories of metabolic traits was estimated using multilevel models. FINDINGS: Total number of participants included was 6537 (12,543 repeated measures). Lower maternal education was associated with more adverse levels of several atherogenic lipids and key metabolic traits among females at age 7 years, but not males. For instance, SII for very small very-low-density lipoprotein (VLDL) concentrations was 0.16SD (95% CI: 0.01, 0.30) among females and -0.02SD (95% CI: -0.16, 0.13) among males. Between 7 and 25 years, inequalities widened among females and emerged among males particularly for VLDL particle concentrations, apolipoprotein-B concentrations, and inflammatory glycoprotein acetyls. For instance, at 25 years, SII for very small VLDL concentrations was 0.36SD (95% CI: 0.20, 0.52) and 0.22SD (95% CI: 0.04, 0.40) among females and males respectively. INTERPRETATION: Prevention of socioeconomic inequalities in cardiovascular disease risk requires a life course approach beginning at the earliest opportunity, especially among females. FUNDING: The UK Medical Research Council and Wellcome (grant ref: 217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). KON is supported by a Health Research Board (HRB) of Ireland Investigator Led Award (ILP-PHR-2022-008). JB, GDS and KT work in a unit funded by the UK MRC (MC_UU_00011/1 and MC UU 00011/3) and the University of Bristol. OR is supported by a UKRI Future Leaders Fellowship (MR/S03532X/1). These funding sources had no role in the design and conduct of this study. This publication is the work of the authors and KON will serve as guarantor for the contents of this paper.

Methodological approaches, challenges, and opportunities in the application of Mendelian randomisation to lifecourse epidemiology: A systematic literature review.

Diseases diagnosed in adulthood may have antecedents throughout (including prenatal) life. Gaining a better understanding of how exposures at different stages in the lifecourse influence health outcomes is key to elucidating the potential benefits of disease prevention strategies. Mendelian randomisation (MR) is increasingly used to estimate causal effects of exposures across the lifecourse on later life outcomes. This systematic literature review explores MR methods used to perform lifecourse investigations and reviews previous work that has utilised MR to elucidate the effects of factors acting at different stages of the lifecourse. We conducted searches in PubMed, Embase, Medline and MedRXiv databases. Thirteen methodological studies were identified. Four studies focused on the impact of time-varying exposures in the interpretation of "standard" MR techniques, five presented methods for repeat measures of the same exposure, and four described methodological approaches to handling multigenerational exposures. A further 127 studies presented the results of an applied research question. Over half of these estimated effects in a single generation and were largely confined to the exploration of questions regarding body composition. The remaining mostly estimated maternal effects. There is a growing body of research focused on the development and application of MR methods to address lifecourse research questions. The underlying assumptions require careful consideration and the interpretation of results rely on select conditions. Whilst we do not advocate for a particular strategy, we encourage practitioners to make informed decisions on how to approach a research question in this field with a solid understanding of the limitations present and how these may be affected by the research question, modelling approach, instrument selection, and data availability.

Trial arm outcome variance difference after dropout as an indicator of missing-not-at-random bias in randomized controlled trials.

Randomized controlled trials (RCTs) are vulnerable to bias from missing data. When outcomes are missing not at random (MNAR), estimates from complete case analysis (CCA) and multiple imputation (MI) may be biased. There is no statistical test for distinguishing between outcomes missing at random (MAR) and MNAR. Current strategies rely on comparing dropout proportions and covariate distributions, and using auxiliary information to assess the likelihood of dropout being associated with the outcome. We propose using the observed variance difference across trial arms as a tool for assessing the risk of dropout being MNAR in RCTs with continuous outcomes. In an RCT, at randomization, the distributions of all covariates should be equal in the populations randomized to the intervention and control arms. Under the assumption of homogeneous treatment effects and homoskedastic outcome errors, the variance of the outcome will also be equal in the two populations over the course of follow-up. We show that under MAR dropout, the observed outcome variances, conditional on the variables included in the model, are equal across trial arms, whereas MNAR dropout may result in unequal variances. Consequently, unequal observed conditional trial arm variances are an indicator of MNAR dropout and possible bias of the estimated treatment effect. Heterogeneous treatment effects or heteroskedastic outcome errors are another potential cause of observing different outcome variances. We show that for longitudinal data, we can isolate the effect of MNAR outcome-dependent dropout by considering the variance difference at baseline in the same set of patients who are observed at final follow-up. We illustrate our method in simulation for CCA and MI, and in applications using individual-level data and summary data.

Association between perinatal mortality and morbidity and customised and non-customised birthweight centiles in Denmark, Finland, Norway, Wales, and England: comparative, population based, record linkage study.

Objectives: To compare the risk of adverse perinatal outcomes according to infants who are born small for gestational age (SGA; <10th centile) or large for gestational age (LGA; >90th centile), as defined by birthweight centiles that are non-customised (ie, standardised by sex and gestational age only) and customised (by sex, gestational age, maternal weight, height, parity, and ethnic group). Design: Comparative, population based, record linkage study with meta-analysis of results. Setting: Denmark, Finland, Norway, Wales, and England (city of Bradford), 1986-2019. Participants: 2 129 782 infants born at term in birth registries. Main outcome measures: Stillbirth, neonatal death, infant death, admission to neonatal intensive care unit, and low Apgar score (<7) at 5 minutes. Results: Relative to those infants born average for gestational age (AGA), both SGA and LGA births were at increased risk of all five outcomes, but observed relative risks were similar irrespective of whether non-customised or customised charts were used. For example, for SGA versus AGA births, when non-customised and customised charts were used, relative risks pooled over countries were 3.60 (95% confidence interval 3.29 to 3.93) versus 3.58 (3.02 to 4.24) for stillbirth, 2.83 (2.18 to 3.67) versus 3.32 (2.05 to 5.36) for neonatal death, 2.82 (2.07 to 3.83) versus 3.17 (2.20 to 4.56) for infant death, 1.66 (1.49 to 1.86) versus 1.54 (1.30 to 1.81) for low Apgar score at 5 minutes, and (based on Bradford data only) 1.97 (1.74 to 2.22) versus 1.94 (1.70 to 2.21) for admission to the neonatal intensive care unit. The estimated sensitivity of combined SGA or LGA births to identify the three mortality outcomes ranged from 31% to 34% for non-customised charts and from 34% to 38% for customised charts, with a specificity of 82% and 80% with non-customised and customised charts, respectively. Conclusions: These results suggest an increased risk of adverse perinatal outcomes of a similar magnitude among SGA or LGA term infants when customised and non-customised centiles are used. Use of customised charts for SGA/LGA births-over and above use of non-customised charts for SGA/LGA births-is unlikely to provide benefits in terms of identifying term births at risk of these outcomes.

Time-sensitive testing pressures and COVID-19 outcomes: are socioeconomic inequalities over the first year of the pandemic explained by selection bias?

BACKGROUND: There are many ways in which selection bias might impact COVID-19 research. Here we focus on selection for receiving a polymerase-chain-reaction (PCR) SARS-CoV-2 test and how known changes to selection pressures over time may bias research into COVID-19 infection. METHODS: Using UK Biobank (N = 420,231; 55% female; mean age = 66.8 [SD = 8·11]) we estimate the association between socio-economic position (SEP) and (i) being tested for SARS-CoV-2 infection versus not being tested (ii) testing positive for SARS-CoV-2 infection versus testing negative and (iii) testing negative for SARS-CoV-2 infection versus not being tested. We construct four distinct time-periods between March 2020 and March 2021, representing distinct periods of testing pressures and lockdown restrictions and specify both time-stratified and combined models for each outcome. We explore potential selection bias by examining associations with positive and negative control exposures. RESULTS: The association between more disadvantaged SEP and receiving a SARS-CoV-2 test attenuated over time. Compared to individuals with a degree, individuals whose highest educational qualification was a GCSE or equivalent had an OR of 1·27 (95% CI: 1·18 to 1·37) in March-May 2020 and 1·13 (95% CI: 1.·10 to 1·16) in January-March 2021. The magnitude of the association between educational attainment and testing positive for SARS-CoV-2 infection increased over the same period. For the equivalent comparison, the OR for testing positive increased from 1·25 (95% CI: 1·04 to 1·47), to 1·69 (95% CI: 1·55 to 1·83). We found little evidence of an association between control exposures, and any considered outcome. CONCLUSIONS: The association between SEP and SARS-CoV-2 testing changed over time, highlighting the potential of time-specific selection pressures to bias analyses of COVID-19. Positive and negative control analyses suggest that changes in the association between SEP and SARS-CoV-2 infection over time likely reflect true increases in socioeconomic inequalities.

Epigenetic clocks and research implications of the lack of data on whom they have been developed: a review of reported and missing sociodemographic characteristics.

Epigenetic clocks are increasingly being used as a tool to assess the impact of a wide variety of phenotypes and exposures on healthy ageing, with a recent focus on social determinants of health. However, little attention has been paid to the sociodemographic characteristics of participants on whom these clocks have been based. Participant characteristics are important because sociodemographic and socioeconomic factors are known to be associated with both DNA methylation variation and healthy ageing. It is also well known that machine learning algorithms have the potential to exacerbate health inequities through the use of unrepresentative samples - prediction models may underperform in social groups that were poorly represented in the training data used to construct the model. To address this gap in the literature, we conducted a review of the sociodemographic characteristics of the participants whose data were used to construct 13 commonly used epigenetic clocks. We found that although some of the epigenetic clocks were created utilizing data provided by individuals from different ages, sexes/genders, and racialized groups, sociodemographic characteristics are generally poorly reported. Reported information is limited by inadequate conceptualization of the social dimensions and exposure implications of gender and racialized inequality, and socioeconomic data are infrequently reported. It is important for future work to ensure clear reporting of tangible data on the sociodemographic and socioeconomic characteristics of all the participants in the study to ensure that other researchers can make informed judgements about the appropriateness of the model for their study population.

Childhood trajectories of internalising and externalising problems associated with a polygenic risk score for neuroticism in a UK birth cohort study.

Background: Neuroticism represents a personality disposition towards experiencing negative emotions more frequently and intensely. Longitudinal studies suggest that neuroticism increases risk of several psychological problems. Improved understanding of how this trait manifests in early life could help inform preventative strategies in those liable to neuroticism. Methods: This study explored how a polygenic risk score for neuroticism (NEU PRS) is expressed from infancy to late childhood across various psychological outcomes using multivariable linear and ordinal regression models. In addition, we employed a three-level mixed-effect model to characterise child internalising and externalising trajectories and estimate how a child PRS associated with both their overall levels and rates of change in 5279 children aged 3-11 in the Avon Longitudinal Study of Parents and Children cohort. Results: We found evidence that the NEU PRS was associated with a more emotionally sensitive temperament in early infancy in addition to higher emotional and behavioural problems and a higher risk of meeting diagnostic criteria for a variety of clinical disorders, particularly anxiety disorders, in childhood. The NEU PRS was associated with overall levels of internalising and externalising trajectories, with a larger magnitude of association on the internalising trajectory. The PRS was also associated with slower rates of reduction of internalising problems across childhood. Conclusions: Our findings using a large, well-characterised birth cohort study suggest that phenotypic manifestations of a PRS for adult neuroticism can be detected as early as in infancy and that this PRS associates with several mental health problems and differences in emotional trajectories across childhood.